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Abstract 


This  report  develops  the  equations  of  motion  for  a  differential  roll  projectile  configuration 
with  seven  degrees  of  freedom.  The  dynamic  equations  are  generated  generically  such  that  the 
forward  and  aft  components  are  mass  unbalanced.  A  hydrodynamic  bearing  exists  between  the 
forward  and  aft  components,  which  couples  the  roll  degree  of  freedom.  A  simulation 
investigation  shows  that  bearing  resistance  and  forward/aft  body  mass  ratio  are  the  dominant 
factors  in  determining  the  roll  dynamics.  For  spin  rates  typical  of  fin-stabilized  projectiles,  the 
trajectory  is  essentially  independent  of  both  bearing  resistance  and  mass  ratio. 
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1.  Introduction 


Compared  to  conventional  munitions,  smart  munitions  involve  more  design  requirements  due 
to  additional  sensors  and  control  mechanisms.  These  additional  components  must  seek  to 
minimize  the  weight  and  space  impact  on  the  overall  projectile  design  so  that  desired  target 
efifects  can  still  be  achieved  with  the  weapon.  The  inherent  design  conflict  between  standard 
projectile  design  considerations  and  new  requirements  imposed  by  sensors  and  control 
mechanisms  has  led  designers  to  consider  more  complex  geometric  configurations.  One  such 
configuration  is  the  differential  roll  projectile.  This  projectile  configuration  is  comprised  of 
forward  and  aft  components.  The  forward  and  aft  components  are  connected  through  a  bearing, 
which  allows  the  forward  and  aft  portions  of  the  projectile  to  spin  at  different  rates.  Figure  1 
shows  a  schematic  of  the  differential  roll  projectile  configuration. 


Typical  flight  mechanic  analysis  of  a  conventional,  single-body,  munition  models  the 
projectile  with  six  degrees  of  freedom.  Dynamic  modeling  of  a  differential  roll  projectile  adds 
an  additional  roll  degree  of  freedom  to  the  equations  of  motion.  This  report  begins  with  the 
development  of  a  dynamic  model  of  a  differential  roll  projectile  in  atmospheric  flight,  including 
the  additional  roll  degree  of  freedom.  The  model  is  derived  such  that  both  the  forward  and  aft 
bodies  can  be  mass  unbalanced.  A  hydrodynamic  bearing  couples  the  forward  and  aft 
components  in  the  roll  axis.  The  mathematical  model  is  utilized  to  show  trends  in  system 
response  as  a  function  of  mass  ratio  and  bearing  resistance. 


2.  Differential  Roll  Projectile  Dynamic  Model 


The  mathematical  model  describing  the  motion  of  the  differential  roll  projectile  allows  for 
three  translation  and  four  rotation  rigid-body  degrees  of  freedom.  The  translation  degrees  of 
freedom  are  the  three  components  of  the  mass  center  position  vector.  The  rotation  degrees  of 
freedom  are  the  Euler  yaw  and  pitch  angles,  as  well  as  the  forward  body  roll  and  aft  body  roll 
angles.  The  equations  presented  here  use  the  ground  surface  as  an  inertial  reference  frame  [1]. 

Development  of  the  kinematic  and  dynamic  equations  of  motion  is  aided  by  the  use  of  an 
intermediate  reference  frame.  The  sequence  of  rotations  from  the  inertial  frame  to  the  forward 
and  aft  bodies  consists  of  a  set  of  body-fixed  rotations  that  are  ordered:  yaw,  pitch,  and 
forward/aft  body  roll.  The  fixed-plane  reference  frame  is  defined  as  the  intermediate  frame 
before  roll  rotation.  The  fixed-plane  frame  is  conunon  to  both  the  forward  and  aft  bodies. 

Equation  (1)  is  the  translation  kinematic  differential  equations  that  relate  time  derivatives  of 
the  mass  center  position  components  to  the  mass  center  velocity  components  in  fixed-plane 
reference  frame. 
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Equation  (2)  is  the  rotation  kinematic  differential  equations  that  relate  time  derivatives  of  the 
Euler  angles  with  angular  velocity  components  in  the  fixed-plane  reference  frame. 
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Equation  (3)  is  the  translation  kinetic  differential  equations  described  in  the  fixed-plane 
reference  fi’ame. 
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Equation  (4)  is  the  rotation  kinetic  differential  equations  described  in  the  fixed-plane 
reference  fi-ame. 
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A  derivation  of  equation  (4),  along  with  definitions  of  the  right  side  components,  is  provided  in 
Appendices  A  and  B. 

As  shown  in  equation  (5),  the  total  applied  force  on  the  complete  configuration  is  provided 
by  the  weight  of  both  the  forward  and  aft  bodies  (w)  and  air  loads  (a). 
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The  weight  portion  of  the  external  loads  is  given  by  equation  (6), 
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while  the  aerodynamic  force  contribution  is  given  by  equation  (7), 
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The  longitudinal  and  lateral  aerodynamic  angles  of  attack  are  computed  using  equation  (8). 


a  =  tan  * 


and  /3  =  tan  * 


rv\ 


(8) 


The  aerodynamic  coefficients  in  equation  (7)  are  functions  of  local  Mach  number  at  the 
projectile  mass  center.  They  are  computed  using  linear  interpolation  from  a  table  of  data.  The 
aerodynamic  forces  and  moments  are  assumed  to  act  solely  on  the  forward  body. 


The  right  side  of  the  rotation  kinetic  equations  contains  the  externally  applied  moments  on 
both  the  forward  and  aft  bodies.  The  external  moment  components  on  the  forward  body  are 
given  by  equation  (9)  and  contain  contributions  from  steady  (sa)  and  unsteady  (ua) 
aerodynamics. 
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The  steady  body  aerodynamic  moment  is  computed  by  a  cross  between  the  distance  vector  from 
the  center  of  gravity  to  the  center  of  pressure,  and  the  steady  body  aerodynamic  force  vector.  The 
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The  unsteady  body  aerodynamic  moment  provides  a  damping  source  for  projectile  angular 
motion  and  is  given  by  equation  (10). 
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Air  density  is  computed  using  the  center  of  gravity  position  of  the  projectile  using  the  standard 
atmosphere  [2]. 

3.  Simulation  Example 


hi  order  to  exercise  the  math  model  discussed  previously,  consider  a  6-ft  long,  120-lb 
projectile.  The  forward  body  is  fin  stabilized  and  the  aft  body  is  an  internal  circular  cylinder. 
Aerodynamic  forces  and  moments  act  on  the  forward  body  only.  For  this  simulation  set,  initial 
forward  body  velocity  is  750  m/s  and  initial  gun  elevation  is  45°.  All  other  states  variables  are 
initially  equal  to  0.  The  projectile  fins  are  canted  slightly  to  provide  a  slowly  rolling  projectile  in 
steady  state. 

Figures  2-1 1  show  the  state  variables  of  the  system  vs.  time  for  the  conditions  mentioned 
above.  The  mass  ratio  of  the  aft  body  to  the  forward  body  is  1%.  Under  these  circumstances,  the 
projectile  has  a  range  of  approximately  18  km.  Cross  range,  yaw  angle,  side  velocity,  vertical 
velocity,  pitch  rate,  yaw  rate,  and  aerodynamic  angle  of  attack  remain  small  throughout  the 
event.  Pitch  attitude  steadily  decreases  from  45°  to  just  below  -60°  at  impact.  Figures  2-11 
remain  the  same,  independent  of  the  bearing  resistance  coefficient.  Figures  12  and  13  show  the 
roll  angle  and  roll  rate  response  as  a  function  of  bearing  resistance  coefficient.  Values  of  bearing 
resistance  coefficient  are  0.000001,  0.000005,  0.00001,  0.00005,  0.0001,  0.0005,  0.001,  0.005, 
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Figure  10.  Yaw  Rate. 
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Figure  11.  Aerodynamic  Angle  of  Attack. 
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Figure  12.  Roll  Angle  (Mass  Ratio  =  1%,  Damping  Coefficient  =  0.01-0.000001). 
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Figure  13.  Roll  Rate  (Mass  Ratio  =  1%,  Damping  Coefficient  =  0.01-0.000001). 
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and  0.01  ft-lbf/rps.  In  Figure  12,  the  lowest  trace  is  the  response  of  the  aft  body  for  the  lowest 
value  of  bearing  resistance.  The  upper  trace  is  the  forward  body  roll  response.  For  a  bearing 
resistance  coefficient  of  0.00005,  the  aft  body  roll  response  is  essentially  the  same  as  the  forward 
body  since  both  bodies  rapidly  couple  in  the  roll  axis.  Figure  13  shows  the  roll  rate  trace  for  this 
simulation  set.  It  is  interesting  that  for  lower  values  of  bearing  resistance,  the  aft  body  roll  rate 
overshoots  the  forward  body  roll  rate  before  settling. 

Figures  14  and  15  show  the  roll  angle  and  roll  rate  response  of  forward  and  aft  bodies  under 
the  same  conditions  as  the  previous  case,  except  the  mass  ratio  is  now  50%  rather  than  1%. 
While  the  basic  character  of  the  roll  response  is  the  same,  the  aft  body  roll  angle,  and  hence  roll 
rate,  build  up  relatively  slowly  due  to  the  increase  in  aft  body  inertia. 
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Figure  14.  Roll  Angle  (Mass  Ratio  =  50%,  Damping  Coefficient  =  0.01-0.000001). 

Figures  16-19  show  the  response  of  the  system  under  the  same  conditions  as  Figures  14  and 
15,  except  the  initial  roll  rate  of  the  aft  body  is  - 100  rad/s.  Like  the  previous  simulation  results, 
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Figure  17.  Cross  Range  (Mass  Ratio  =  50%,  Damping  Coefficient  =  0.01-0.000001). 


Time  (s) 


Figure  18.  Angle  of  Attack  (Mass  Ratio  =  50%,  Damping  Coefficient  =  0.01-0.000001). 
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Figure  19.  Cross  Range  (Damping  Coefficient  =  0.0005,  Mass  Ratio  =  l%-50%). 

lower  values  of  bearing  resistance  produce  slower  roll  response  in  the  aft  body.  For  larger  splits 
in  the  forward  and  aft  body  roll  rates,  the  trajectory  begins  to  change  as  a  function  of  bearing 
resistance  owning  to  the  fact  that  the  roll  response  is  sensitive  to  bearing  resistance.  In 
particular,  Figure  17  shows  the  cross  range  under  these  circumstances.  While  the  spray  m  the 
trajectory  is  only  on  the  order  of  10  m,  it  points  to  the  fact  that  if  the  forward  and  aft  bodies 
possess  substantially  different  initial  roll  rates,  the  trajectory  becomes  a  function  of  bearing 
resistance. 

Figures  19-22  show  system  response  under  the  same  conditions  as  Figures  16-19,  except 
now  the  mass  ratio  is  varied.  Figure  22  shows  the  roll  rate  response.  Due  to  aft  body  inertia 
changes,  the  roll  response  varies  significantly  with  mass  ratio.  Subsequently,  the  trajectory 
begins  to  vary  as  well.  Similar  to  the  previous  case,  the  trajectory  spray  is  on  the  order  of  10  m; 
this  shows  that  trajectory  of  configurations  with  forward  and  aft  bodies  operating  at  significantly 
different  roll  rates  is  a  function  of  the  mass  ratio. 
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Figure  20.  Roll  Angle  (Damping  Coefficient  =  0.0005,  Mass  Ratio  =  l%-50%). 
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Figure  21.  Side  Velocity  (Damping  Coefficient  =  0.0005,  Mass  Ratio  =  l%-50%) 
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Figure  22.  Roll  Rate  (Damping  Coefficient  =  0.0005,  Mass  Ratio  =  l%-50%). 


4.  Conclusions 


The  equations  of  motion  for  a  differential  roll  projectile  configuration  with  seven  degrees  of 
freedom  have  been  developed  and  exercised.  The  dynamic  equations  allow  the  forward  and  aft 
bodies  to  be  mass  unbalanced.  A  hydrodynamic  bearing  between  the  forward  and  aft 
components  couples  the  roll  degrees  of  freedom.  Bearing  resistance  and  forward/aft  body  mass 
ratio  are  the  dominant  factors  in  determining  the  roll  dynamics.  For  spin  rates  typical  of 
fin-stabilized  projectiles,  the  trajectory  is  essentially  independent  of  both  bearing  resistance  and 
mass  ratio.  However,  for  configurations  with  the  forward  and  aft  components  operating  at 
significantly  different  roll  rates,  the  trajectory  depends  on  the  mass  ratio  and  bearing  resistance. 


17 


Intentionally  left  blank. 


18 


5.  References 


1.  Etkin,  B.  Dynamics  of  Atmospheric  Flight.  New  York;  John  Wiley  and  Sons,  1972. 

2.  VonMises,  R.  Theory  of  Flight.  New  York:  Dover  Publications  Inc.,  1959. 

3.  Close,  C.  M.,  and  D.  K.  Frederick.  Modeling  and  Analysis  of  Dynamic  Systems.  New  York: 
John  Wiley  and  Sons,  1995. 


19 


Intentionally  left  blank. 


20 


Appendix  A: 

Constraint  Forces  and  Moments 
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The  rotation  kinetic  differential  equations  are  derived  by  first  splitting  the  two  body  system  at 
the  bearing  connection  point  Figures  A-1  and  A-2  show  the  external  loads  and  internal 
constraint  forces  acting  on  both  the  forward  and  aft  bodies,  respectively. 
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Figure  A-1.  Forces  and  Moments  on  the  Forward  Body. 
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Figure  A-2.  Forces  and  Moments  on  the  Aft  Body. 

The  constraint  force,  ,  and  the  constraint  moment  M^. ,  couple  the  forward  and  aft  bodies. 

Key  to  the  development  of  the  rotation  kinetic  differential  equations  is  the  ability  to  solve  for  the 
constraint  forces  and  moments  as  a  function  of  state  variables  and  time  derivatives  of  state 
variables. 
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An  expression  for  the  constraint  force  can  be  obtained  by  subtracting  the  translation  dynamic 
equations  for  both  bodies. 
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The  acceleration  of  the  mass  center  of  the  forward  and  aft  bodies,  5^/,  and  a^u ,  can  be 
expressed  in  terms  of  the  acceleration  of  the  composite  body  mass  center.  After  making  this 
substitution,  the  constraint  force  components  in  the  fixed-plane  reference  frame  can  be  expressed 
in  the  following  marmer: 
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The  matrices  [Fp],  [F^],  and  {F^}  are  complicated  functions  of  the  state  variables  and  the 
geometry  of  the  configuration. 

The  components  in  the  fixed-plane  reference  frame  of  the  moment  of  the  constraint  force, 
acting  on  the  forward  body  about  the  forward  body  mass  center,  can  be  written  in  the  following 
manner 
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In  a  similar  way,  the  components  in  the  fixed-plane  reference  frame  of  the  moment  of  the 
constraint  force,  acting  on  the  aft  body  about  the  aft  body  mass  center,  can  be  written  in  the 
following  manner: 
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The  matrices  [M^],  [M^],  and  are  also  complicated 

functions  of  the  state  variables  and  the  geometry  of  the  configuration. 
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Appendix  B: 

Rotation  Kinetic  Equations 
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The  rotation  kinetic  differential  equations  are  derived  by  first  writing  the  Euler  equations  for 
the  forward  and  aft  bodies  separately.  These  equations  are  expressed  in  the  fixed-plane  reference 
frame  are  general,  and  allow  for  a  fully  populated  inertia  matrix  and  mass  unbalance.  Equations 
(A-2),  (A-3),  and  (A-4)  are  substituted  into  both  sets  of  rotation  kinetic  equations  for  the  forward 
and  aft  bodies.  At  this  point,  both  sets  of  equations  still  have  unknown  constraint  moments  at  the 
bearing  connection  point.  To  eliminate  the  bearing  constraint  moments  in  the  y  and  z  direction 
in  the  fixed-plane  coordinate  system,  the  y  and  z  components  of  the  rotation  kinetic  equations  for 
the  forward  and  aft  bodies  are  added  together  to  form  two  dynamic  equations  that  are  free  of 
constraint  moments.  In  this  way,  the  constraint  moments  at  the  bearing  have  been  eliminated 
analytically. 

The  forward  and  aft  bodies  are  connected  through  a  hydrodynamic  bearing.  The  moment 
transmitted  across  a  hydrodynamic  bearing  can  be  modeled  as  viscous  damping.^  The 
constitutive  relation  governing  the  constraint  moment  transmitted  across  a  hydrodynamic  bearing 
is  given  by  Equation  (B-1). 


My  =Cy(pp-p^).  (B-1) 

If  the  viscous  damping  coefficient,  ,  equals  zero,  then  the  forward  and  aft  body  connection  is 
frictionless. 

The  effective  inertia  matrix  is  a  4  x  4  matrix  that  is  a  combination  of  the  inertia  matrices  of 
both  the  forward  and  aft  bodies.  As  an  aid  in  developing  a  formula  for  the  effective  inertia 
matrix,  define  the  following  intermediate  matrices: 

(B-2) 


^  Close,  C.  M.,  and  D.  K.  Frederick.  Modeling  and  Analysis  of  Dynamic  Systems.  New  York;  John  Wiley  and 
Sons,  1995. 
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Using  Equations  (B-2),  (B-3),  (B-4),  and  (B-5),  elements  of  the  effective  inertia  matrix  can  now 
be  formed. 
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•^3.3  ~  ■^f7-2,2  ■^^.2  ■^^■^2.2  ''’  ^ AF2,i  ’ 


(B-16) 


■^3.4  “  ■•■  ^Al2.3  '*’  ^^2,3 


(B-17) 


■^4,1  “  ■*■  » 


(B-18) 


•^4,2  ~  -^FAs.,  » 


(B-19) 


■^4,3  ~  ^FF3.2  ^AA3.2  ■*■  -^FAj.!  ■*■  ^AF3.2  ’ 


(B-20) 


^4,4  -  -^FFj  j  ■*■  ^AA3.3  +  ^FA3.3  -^AFs,,  * 


(B-21) 


The  elements  of  the  right-hand  side  vector  are  given  by  Equations  (B-22)  and  (B-23). 


gi^  [Mpx  I  (Pf 

gp^  ►  =  •  MpY  ’“ISfJ  Q 
SFj  .^fz.  .  ^  . 


(B-22) 


gA.  [Max'!  [PaI 

gA,  f  =  lMAYr~PA|  q  K'^Ao}- 


[gA3j  [^AZ 


(B-23) 


The  matrices  [5^]  and  [5^]  in  Equations  (B-22)  and  (B-23)  are  given  by  Equations  (B-24)  and 


(B-25)  as  follows: 


31 


[Sf  ]=  [7>  Y  [h  I’"  [Sf  ]  [If  ][7>  I  and 


(B-24) 


[sj=  [r.F  [/.f [rjqs  J[/ J[rJ, 


where 


[S,]  = 


0  s.  q-c.  r  c.  q  +  s.  r 

VF  ^  9f  VF^  VF 

%r-s^^q  0  -pp 


■s^r 


Pf 


0 


[5J= 


0 

0  -p^ 

Pa  0 


[7>]  = 


0  0 


"Of  ''Of 


[7’J  = 


1 

0  -5 


0  0 


Oa 


"Oa 


'Oa. 


[t,]  = 


0 


0 


<PF+Ie^)s^,  (Pf+Ib^K, 
■(PF+Ie’'KF  -(PF+Ie’')s^, 


,  and 


[TA  = 


0  0  0 
0  -(PA  +  ter)s^,  (PA+l0r)c^^ 

0  -  (Pa  +  Ier)c^,  “  (Pa  + 


(B-25) 


(B-26) 


(B-27) 


(B-28) 


(B-29) 


(B-30) 


(B-31) 
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List  of  Symbols 


x,y,z 


e,y/ 

<t>F 

u,v,w 


Pf 


Pa 


X,Y,Z 

Lp,Mp,Np 

nip 

niA 

m 


[/a] 

[/] 

D 

C, 

a 

P 

[Tp] 

[TJ 

Cy 

V 


Position  vector  components  of  the  center  of  mass  expressed  in  the  inertial 

reference  frame 

Euler  pitch,  yaw  angles 

Euler  roll  angle  of  the  forward  body 

Euler  roll  angle  of  the  aft  body 

Translation  velocity  components  of  the  center  of  mass  resolved  in  the 
fixed-plane  reference  frame 

Roll  axis  component  of  the  angular  velocity  vector  of  the  forward  body 
expressed  in  the  fixed-plane  reference  frame 

Roll  axis  component  of  the  angular  velocity  vector  of  the  aft  body  expressed  in 
the  fixed-plane  reference  frame 

Components  of  the  angular  velocity  vector  of  both  the  forward  and  aft  bodies 
expressed  in  the  fixed-plane  reference  frame 

Total  external  force  components  on  the  projectile  expressed  in  the  fixed-plane 
reference  frame 

External  moments  on  the  forward  body  expressed  in  the  fixed-plane  reference 
frame 

External  moments  on  the  aft  body  expressed  in  the  fixed  plane  reference  frame 
Forward  body  mass 
Aft  body  mass 
Total  projectile  mass 

Mass  moment  of  inertia  matrix  of  the  forward  body  with  respect  to  the  forward 
body  reference  frame 

Mass  moment  of  inertia  matrix  of  the  aft  body  with  respect  to  the  aft  body 

reference  frame 
Effective  inertia  matrix 
Projectile  characteristic  length 
Projectile  aerodynamic  coefficients 

Dynamic  pressure  at  the  projectile  mass  center 

Longitudinal  aerodynamic  angle  of  attack 
Lateral  aerodynamic  angle  of  attack 

Transformation  matrix  from  the  fixed-plane  reference  frame  to  the  forward  body 
reference  frame 

Transformation  matrix  from  the  fixed-plane  reference  frame  to  the  aft  body 

reference  frame 

Viscous  damping  coefficient 

Magnitude-of-mass  center  velocity 
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